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An experimental investigation has been made of the effects of wall perturbations in 
multimode rectangular waveguides, whose dimensions were such that the modes TE no for 
71=1,2 ... 15 were propagating modes and all of the others were cutoff. Measurements 
were made at X-band frequencies. The investigation has shown that when the lowest 
ordered mode is incident upon a perturbation in the waveguide wall, all of the higher ordered 
propagating modes will be excited in both directions in the waveguide, but the distribution 
of energy among the modes is critically dependent upon the shape and size of the perturba- 
tion. 

Measurements on this waveguide with various wall perturbations may be used to gain 
insight into the effects of ionospheric perturbations on terrestrial VLF propagation. 



1. Introduction 

Perturbations in the walls of waveguides can in 
the general case excite all orders of modes traveling 
away from the perturbation. In the case of a wave- 
guide propagating only a single mode, the higher 
modes are of little importance except in the immedi- 
ate vicinity of the perturbation. However, in a 
waveguide which propagates more than one mode, 
several of the higher modes excited by a perturbation 
can be significant at great distances from it. in 
such a waveguide the relative amplitudes and phases 
of the various modes are dependent upon the shape 
and size of the perturbation. Because of this de- 
pendence it is possible, in theory at least, to gain 
information about the size and shape of a perturba- 
tion from measurements of the various propagating 
modes. This principle might therefore be used to 
gain information about perturbations of the iono- 
sphere from measurements of VLF propagation in 
the earth-ionosphere waveguide. It is felt that some 
insight into the effects of ionospheric perturbations 
upon VLF propagation can be acquired by the 
laboratory study of the effects of wall perturbations 
in a microwave waveguide having a width, in wave- 
lengths at the microwave frequency of operation, 
equal to the height of the earth-ionosphere wave- 
guide in wavelengths at the VLF frequency of 
interest. The boundaries of the microwave wave- 
guide should be chosen so that their reflection co- 
efficients are the same as those in the earth-ionosphere 
waveguide. It is theoretically possible to select the 
materials and dimensions of the microwave wave- 
guide so that it is an exact scale model of the earth- 
ionosphere waveguide; then the electromagnetic 
behavior of the model will precisely give the electro- 
magnetic behavior of the earth-ionosphere wave- 
guide. The technical problems involved in the 
construction of an exact scale model are formidable; 
therefore it was decided that the initial investiga- 
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tions would be made with a simple waveguide of 
rectangular cross section. Such a waveguide bears 
little resemblance to the earth-ionosphere waveguide; 
it transmits power along an axis rather than between 
two spheres as in the earth-ionosphere waveguide. 
However, if its width is properly chosen and the 
proper wall materials are chosen, then its behavior 
for certain wall perturbations should be similar to 
that of the earth-ionosphere waveguide. J. R. Wait 
has recently given strong theoretical justification for 
this sort of experimental study [Wait, 1903]. A 
further important use of such a waveguide is that of 
checking the accuracy of various approximate 4 theo- 
retical solutions for fields in waveguide having wall 
perturbations [Wait and Murphy, 1957; Wait, 1960a, 
1960b, 1961, 1962; Morgan, 1950]. 



2. Theory 

The earth can usually be considered as a good 
conductor at VLF, therefore it is a straightforward 
matter to represent it in an electromagnetic model. 
It is simply necessary to use a material having a 
conductivity equal to that of the earth multiplied 
by an appropriate scale factor. Furthermore, since 
the conductivity is scaled upward when the size 
is scaled downward it follows that, in a model small 
enough to be placed in a laboratory, the conductivity 
of the material representing the earth would be very 
high, so high, in fact, that a perfect conductor can 
be used with very little error. 

On the other hand, the representation of the 
onosphere in the model is somewhat more difficult. 
The magnitude of the refractive index of the iono- 
sphere at VLF under certain rather common condi- 
tions is approximately unity; it can range somewhat 
above or below this figure. It is more difficult to 
model a medium with a refractive index <0 than 
one with a refractive index >1. Therefore it was 
decided to first model the case for which the refractive 
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index, A 7 , is defined by 2 

where L is a positive number ranging from about 
0.1 to 3. The surface impedance, Z, corresponding 
to this refractive index is 



«-^(-"-0' 



where 



L 

7?o^l207r; 



and 6 is the angle of incidence. For the lower 



ordered modes, the ansie of incidence is close to 



2' 



and, for such angles the surface impedance, as given 
above, has an imaginary part which is negative for 
large L and positive for small L. A suitable model, 
therefore, would be as shown in figure 1. The 
bottom wall would be a very good conductor, and 
the top wall would be a material having a surface 
impedance given by the above equation for Z. The 
desired modes of operation of the waveguide in 
figure 1 are TM modes which have no tangential 
components of the magnetic field intensity, H, on 
the side walls and no variation of the field quantities 
in the direction of the normal to the broad walls. 
This would require that the broad walls have an 
essentially infinite permeability. Such walls are 
difficult to even approximate; therefore it was 
decided that the type of wave guide shown in figure 1 
is unsuitable. J. R. Wait has suggested a way to 
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MATERIAL TO SIMULATE 
THE IONOSPHERE 



IMPERFECT CONDUCTOR 
TO SIMULATE THE GROUND 



Figure 1. Cross section of a possible laboratory waveguide, 
a two-dimensional model of the earth-ionosphere region. 

2 The rationalized M.K.S. systemof units will be used. The time dependence 
is assumed to be e to,i , where i= V— 1, co is angular frequency and t is time. 



overcome the disadvantages of the waveguide of 
figure 1 [Maley and Bahar, 1963]. He suggested 
the one shown in figure 2. The part of the waveguide 
above the plane of symmetry in figure 2 may be 
considered as an electromagnetic dual of the wave- 
guide in figure 1 for the case in which the bottom 
wall in figure 1 is a perfect conductor. The term 
electromagnetic dual is used to mean that the 
fields and the parameters are interchanged according 
to the relations 



where, in general, e and p are complex. Maxwell's 
equations are invariant to this transformation; so 
fields which exist in an electromagnetic system can 
be directly related to those which exist in the dual 
of the system. Thus if the field is measured in the 
waveguide of figure 2, the field in the waveguide of 
figure 1 can be calculated from the measurements by 
the above transformation. Using this transforma- 
tion, the bottom wall in figure 1 (this wall is assumed 
to be a perfect conductor) transforms into a plane 
along which there can be no tangential component 
of H. But this is exactly the behavior of the plane 
of symmetry of figure 2 when only TE n0 modes 
having E symmetric about the plane of symmetry 
are considered. In a completely symmetric system 
only these modes exist, and therefore the plane of 
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Figure 2. Cross section of the laboratory waveguide used in 
the experimental investigation. 
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symmetry may be regarded as the exact dual of a 
perfect conductor. The top wall in figure 1 trans- 
forms into a wall having a surface admittance, Y, 
equal to 

V( w (cos 2 0-^Z) 1/2 
r) (l — lL) 

The broad walls of figure 1 are transformed into 
perfect conductors. The waveguide in figure 2 can 
be approximated to a good order of accuracy using 
readily available materials. The broad walls can 
be made of any metallic substance and the top and 
bottom walls can be made of commercially available 
imperfectly conducting materials. Since it is desired 
that there be no variations of the field quantities 
in the horizontal direction in figure 2, it is advisable 
to make the narrow dimension of the waveguide 
sufficiently small that only TE, l0 modes can propa- 
gate. For the lower ordered modes in the earth 
ionosphere waveguide, the reflection coefficient of t he 
ionosphere, R u given by [Wait, 1900a] 



Ri= 



cos 6—Zt]o 1 

COS 6-{-Zr)o l 



can be roughly approximated by —1 which in (he 
dual model corresponds to +1 or a perfect conductor. 
Therefore the first investigations were made with 
brass top and bottom walls on the waveguide. 
Later an imperfectly conducting material was used 
to give a better approximation; a still closer ap- 
proximation is now being planned. 

3. Experimental Equipment 

The model waveguide was constructed from three 
8-ft long sections with uniform cross section 25.4 
cmX0.95 cm. The operating frequency was 9-1 9 c/s, 
hence the length of the waveguide was about 225 
A (X being the free space wavelength) representing 
4,500 km in the earth-ionosphere waveguide at 15 
kc/s. The width of the waveguide was slightly over 
7.5 A and the height was less than half a wavelength. 
Symmetry was maintained in the waveguide at all 
times (even when perturbations were inserted) to 
maintain the dual boundary conditions, hence only 
the first eight TE n , modes with symmetric electric 
fields (having n only odd integers) were propagating 
modes, the others being evanescent. A nonreflecting 
load was used to terminate the waveguide. The 
"perfectly conducting' ' narrow walls of the wave- 
guide represent the idealized ionospheric boundary 
in the earth-ionosphere waveguide (having a vanish- 
ing tangential magnetic field). To simulate the 
physical characteristics of the ionosphere more 
closely, these walls will be padded with a lossy 
dielectric material that is now in the process of being 
manufactured. 

At large distances from a transmitting vertical 
dipole in the earth-ionosphere waveguide the field 
would be represented essentially by a TM 10 mode, 
being the least attenuated of all the propagating 
modes. Hence it was desired to launch a TE 1>0 



mode into the model waveguide. This was obtained 
by transmitting a TE 1>0 mode from a standard 
X-band waveguide through a 2.5 m long taper with 
a sinusoidal profile into the model waveguide. The 
wave thus launched into the model waveguide had 
a higher mode content (essentially TE 3i ) of 8 percent 
approximately. 

Amplitude and phase measurements were carried 
out by probing the electric field along the center 
line of the waveguide and across transverse sections 
of it, through a series of holes, rather than slots. In 
order to observe the detail interference structure 
of the field pattern it was necessary to space the 
holes such that at least eight readings would be 
taken per free space wavelength along the length of 
the waveguide. It was necessary to take a minimum 
of 500 readings on either side of the pert urbed region 
to recognize the presence of all the higher modes 
propagating in the waveguide. 

A block diagram of the amplitude and phase 
measuring system is given in figure 3. Since exten- 
sive amplitude data was to be obtained, it was 
desired to have this quantity plotted directly. The 
amplitude modulated signal sampled by a probe 
with a square law detector was amplified and 
rectified. This rectified signal was passed through 
a square root circuit before being recorded on an 
X-Y recorder. The X axis representing the posi- 
tion of the probe and the )' axis the amplitude of 
the electric field at the sampling holes. Phase read- 
ings were taken by comparing the phase of the 
sampled electric field with that of a reference signal 
in a Magic Tee with the help of a precision phase 
shifter. A general view of the system with only 
one section connected to the taper section is shown 
in figure 4. 

4. Electric Field Measurements 

The first series of experiments were conducted 
with the uniform cross section of the waveguide 
locally perturbed by sinusoidal shaped slabs (of 
the same material as the waveguide walls) placed 
symmetrically along the walls. The first pair of 
perturbations were approximately 2K X long such 
that the half width of the waveguide varied as 



h(z)=h 



[>- 



a cos 



w(z — z q ) 
2z 



] 
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h is the half width of the uniform waveguide and 
(ha) is the maximum height of the perturbation and 
2z is the length of the perturbation (see fig. 5), 
a=0.2 representing a 20 percent local depression 
in the height of the ionosphere. 

Measurements of the transverse and longitudinal 
variations of the electric field (amplitude and phase) 
were taken on both sides of the perturbation. The 
amplitude measurement along the axis of the guide 
revealed distinctly different standing wave patterns 
on either side of the perturbation (fig. 6a). These 
patterns could be precisely explained theoretically 
in terms of interaction between the principal mode 
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Figure 3. Block diagram of the system used for measurement of the electric field intensity and phase. 
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Figure 4. General view of a section of the model waveguide Figure 5. Axial section of waveguide with sinusoidal 

backed by the amplitude and phase measuring device. perturbation inserted. 
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Figure 6, a, b. Axial amplitude and phase variations with 
sinvsoidal perturbation & = 0.2. 



and the higher order modes generated by the per- 
turbation. While the phase variation along the 
axis of the waveguide was very close to linear, the 
phase anomaly A(p (difference between the phase 
variation along the axis of the waveguide in the 
perturbed case and in the unperturbed case) re- 
markably resembled the amplitude variations (fig. 
6b). This too can be accounted for in terms of mode 
interaction. The transverse amplitude variation 
(fig. 7a) emphasized the higher mode content. The 
very near symmetry of the measured variations 
implied that only the odd modes were being excited 
as is required by the duality relationship between 
the model and the actual problem. The transverse 
phase variations (fig. 7b) also clearly showed up the 
higher mode content. 

To explain the close correspondence between the 
axial amplitude and phase anomaly variations in the 
perturbed waveguide, one notes, for instance, that 




Figure 7, a, b. Transverse amplitude and phase variations 
with sinusoidal perturbation sl = 0.2. 
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Figure 8a. Total electric field represented by higher modes 
added to principal mode. 




Figure 8b. Relationship between amplitude and phase 
anomaly. 



in the region between the generator and the per- 
turbation the resultant field comprises the (forward 
propagating) principal mode and a sum of reflected 
(backward propagating) modes. 

In the schematic representation of the resultant 
field variations along the axis of the waveguide 
(fig. 8a) E x and <p x are the unperturbed amplitude 
and phase of the electric field (principal mode). 
E t and <p are the resultant (perturbed) amplitude 
and phase of the electric field, while E A and A<p are 
the amplitude and phase anomaly. Note that 
while Ei rotates counterclockwise, E& rotates clock- 
wise (towards generator). Hence, while the stand- 
ing wave patterns of E t and A<p have the same 
wavelength, and the variation of E t lags that of 
A<p by 7r/2 radians (fig. 8b). 

While the standing wave pattern in a one mode 
system (in which only one mode propagating in each 
direction is present) has a standing wave wave- 
length of }i \ g , those measured in the perturbed 
(multimode) waveguide were 0.85 X^ in the region 
between she generator and the perturbation and 
about 1.2 \ g in the region beyond the perturbation 
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The total field in the region between the genera- 
tor and the perturbation is approximately 



+B I u e 



\Xgi x»i5/ . r 



71=1,3,5 ... 15, — (3— s )>Xo; 
and in the region beyond the perturbation 

+ . . . fl?e«* GJr£) 2 + . . . 



+Rlle i2 * (53-xs) *} 



n= 



1,3,5,7 . 



15, 



(2-2 )>X , 



where E\ is the incident principal (forward propa- 
gating) mode in region I and E\R I 1 I =El I is the 
principal mode transmitted into region II. R\ 
= E I JE 1 1 and R 1 *=E I n I /E 1 l are the reflected and 
transmitted higher mode conversion coefficients 
respectively, 2 is the distance along the axis of the 
waveguide and \ gn the waveguide wavelength of the 
nth mode. 

The interaction between the 15th mode (back- 
ward propagating) and the principal mode in region 
1 has an effective standing-wave wavelength, Xe ff i 5 ; 
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NowX g „=X [l-(^) 2 ] I/2 =^5 where the free 

space wavelength X =3.33 cm and the half width 
of the waveguide h=12.7 cm. Thus 



X e ff 



COS 61 



COS ^ + cos n 



=0.85 X^. 



Similarly, in region I beyond the perturbation, the 
interaction between the 15th (forward propagating) 
mode and the principal mode has an effective 
standing-wave wavelength 



XfliX, 



g!5 



S?l n n~ — 1.22 X„i* 

8 cos d l — cos 15 



From the above results it seems that for the par- 
ticular type of perturbation introduced into the 
waveguide in the above experiment i? i5 was the 



largest of all the mode conversion coefficients, 
since a predominant Xe ff i 5 = 0.85 \ gl and Xeffi 5 =1.22 \ gi 
were actually obtained from the amplitude and 
phase data taken. Looking into the transverse 
field distributions of the different modes one could 
perhaps draw some qualitative explanation for the 
observation that R 15 is the largest of all the mode 
conversion coefficients for the particular form of 
perturbation used in the above experiment. 

Figure 9 is the pattern of the electric field distri- 
bution across the waveguide for the TE 15 , mode 
in the uniform cross section (0.95 cm X 25.4 cm) 
waveguide. The first null in the electric field is 

— '■=—, '—=1.7 cm from the uniform waveguide wall. 
n 15 

In general, for the nth mode the first null will be 

— — cm from the waveguide wall. 
n 

The average height of the sinusoidal perturbation 

2 
used in the above experiments is -X2.54 cm=1.62 

7T 

cm. Now the transverse electric field pattern of 
the 15th mode in the uniform 25.4 cm waveguide 
could also meet the required boundary condition 
in a 2(12.7 — 1.7) cm uniform waveguide, and in the 
region of the perturbation the average width of the 
guide is 2(12.7—1.62) cm. Hence, this could 
imply that of all the propagating modes in the uni- 
form 25.4 cm waveguide the 15th mode meets the 
boundary condition in the region of the perturba- 
tion most closely. Hence, the 15th mode would 
be the most favored of all the higher modes in the 
region of the perturbation and R l5 the largest of all 
the mode conversion coefficients. 

To support the above suggestion concerning the 
reason for the particular higher mode content 
observed, the experiment was repeated with another 
sinusoidal perturbation about 3X long and with 
(ha) =3.05 cm instead of 2.54 cm as in the previous 
case. In this case the average height of the pertur- 
bation was 1/13 the width of the unperturbed wave- 
guide. The amplitude variations along the axis of 
the guide showed that the TE 13 mode was the most 
predominant mode excited by the perturbation (fig. 
10). The transverse amplitude variations plotted in 
figure 11 also very clearly show the interference of 
the TE 13> mode. 
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Figure 9. Transverse variations of the electric field for the 
TEo,i5 mode showing situation when first null coincides with 
average height of perturbation. 
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Figure 10. Axial amplitude variations with sinusoidal 
perturbation 2l = 0.24- 




Figure 11. 
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Further experiments with sinusoidal perturbations 
showed that for perturbation lengths over 20A ffl , t he- 
mode conversion tends to become negligible if its 
height does not exceed 20 percent of the waveguide 
half width. 

By narrowing a length of the waveguide to a 
uniform 20.3 cm X 0.95 cm cross section, it was 
possible to simulate conditions similar to those 
encountered when the path of propagation is crossed 
by sunrise or sunset. In the case of a night to day 
path it was noted that more of the higher order modes 
were excited at the transition, being both reflected 
from it and transmitted beyond it. For the day to 
night path the lower order modes were the more 
predominant, but it was remarkable to note that in 
this case no matter what the shape of the transition 
joining the narrow waveguide to the broad one all 
the reflections from the transition region in the 
backward direction were negligible while the wave 
pattern beyond the transition region was rich in 
higher mode content as seen in figure 12 and figure 13 
showing the longitudinal and transverse amplitude 
variations. 

The above observations have suggested an approxi- 
mate theoretical approach to the solution of prop- 




DIRECTION OF PROPAGATION 



Figure 12. Axial amplitude variations for day to night 
transition path. 




Figure 13. Transverse amplitude variation* for day to night 

transition path. 



agation problems in in nltimode waveguides yielding 
results that compare favorably with the experimental 
results. These solutions which will be dealt with in 
det ail in a subsequent paper, predict also that the 
mode which is closest to being cutoff in the greater 
part of the transition region of the waveguide, tends 
to be predominant. 

The mode closest to cutoff is that mode for which 



cos 6 n 



(> 



<mi 



is a minimum in the transition region. In the above 
equation h=l2.7 cm has been replaced by h(z) the 
half width of the waveguide as a function of distance- 
in the transition region. These theoretical results 
could also be applied to explaining the mode content 
in the cases when sinusoidal perturbations were 
introduced. 

5. Conclusions 

The experimental investigations thus far have 
shown that perturbations in the walls of multimode 
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waveguides excite all orders of propagating modes 
traveling away from the perturbation. In all cases 
the predominant mode on both sides of the pertur- 
bation has been the lowest ordered mode. This is 
to be expected since energy incident upon the per- 
turbation was predominantly in this mode. In most 
cases there was also a substantial amount of the 3d 
mode (this would correspond to the second mode in 
the earth-ionosphere waveguide 3 ); this was partly 
due to the fact that the incident wave contained 
some 3d mode energy and partly because of conver- 
sion of energy from the lowest mode at the perturba- 
tion. The higher ordered modes observed were 
primarily caused by the presence of the perturbation. 
For the perturbations used it was found that the 
energy in modes 13 or 15 or both was substantially 
greater than that in modes 7, 9, and 11. It was also 
observed that small changes in the width of the per- 
turbation made a large change in the relative ampli- 
tudes of modes 13 and 15. The general conclusion 
is that the distribution of energy among the various 
propagating modes is critically dependent upon the 
size of the perturbation. This is particularly true 
in the case of perturbations less than about three 
wavelengths long. For longer perturbations the 
mode conversion becomes less until at a length of 
20 wavelengths (and a height less than a wavelength) 
there is hardly any noticeable mode conversion at all. 



3 The modes 1, 3, 5, . . . 15 in the model waveguide correspond to modes 
1, 2, 3, . . . 8, respectively, in the earth ionosphere waveguide. 



The measurement techniques have now been re- 
fined to the extent that quantitative comparisons 
can be made with theoretical solutions. It is ex- 
pected that data on such comparison will become 
available soon. 
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